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We report a theoretical study of an overdamped Brownian particle dynamics in the presence of 
both a spatially modulated one-dimensional periodic potential U(x) and a periodic alternating force 
(AF). As the periodic potential U(x) has a low symmetry (a ratchet potential) the Brownian particle 
displays a broken symmetry motion with a nonzero time average velocity. By making use of the 
Green function method and a mapping to the theory of Brillouin bands the probability distribution 
P(x,t) of the particle coordinate x is derived and the nonlinear dependence of the macroscopic 
velocity on the frequency uj and the amplitude n of AF is found. In particular, our theory allows 
to go beyond the adiabatic limit (to = 0) and to explain the peculiar reversal of the velocity sign 
found previously in the numerical analysis. 
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The model of a particle moving in a periodic one- 
dimensional potential U(x) has been used to explain di- 
verse fascinating phenomena in various fields of physics, 
chemistry and biologyEFll Well known examples of 
such phenomena are Brillouin bands and Bloch oscil- 
lations in thev-,solid state theoryu, dynamics of Joseph- 
son junjctionstJ, energy transport in various biological 
systemscl, etc. 

The peculiar property of a particle motion in the pe- 
riodic potential is the presence of two different states, 
a static state and a whirling (dynamic) state. In the 
dynamic state, the particle has a non zero value of the 
time averaged velocity v =< x(t) >. As the system 
is subject to an external noise, the " 'Brownian" parti- 
cle motion becomes more complex (it displays both ran- 
dom damped oscillations and random jumps between the 
potential wellso) , the difference between libration and ro- 
tation states disappears. 

The theory of Brownian particle motion in a specific 
periodic potential of cos-like shape and in the presence 
of a jSonstant driving force (DF) has been developed in 
Refsfl™. As a particular result it was found that in the 
absence of DF the mean particle velocity is zero. More- 
over, it is well knownliiffl that equilibrium noise only can 
not lead to the fluctuation induced transport, and corre- 
spondingly, v — for the particle motion in an arbitrary 
periodic potential U (x) . In this case the diffusive motion 
of particle occurs. The probabilities of fluctuation in- 
duced particle escape from the potential well to the left 
and to the right are identical. These probabilities are de- 
termined by the amplitude of potential U(x) and do not 
depend on the symmetry of potential. 

The situation changes drastically when a Brownian 
particle is subject to both a periodic spatially modulated 
potential (PP) and a periodic alternating force (AF). It 
was shown in Ref. |ll| by numerical analysis of the cor- 
responding Fokker-Planck equation that in this case the 
reflection symmetry of an overdamped particle motion 
(x to —x) can be broken, and directed transport occurs. 



The specific condition of such transport is the low sym- 
metry of the periodic potential U(x) (a ratchet potential). 
Moreover, by making use of the symmetry arguments and 
the numerical analysis it was shown that directed trans- 
port occurs even in the presence of symmetric-PP but as 
the periodic AF has a low symmetry in timea. 

An analytical description of a Brownian particle mo- 
tion in an anisotropic periodic potential has been carried 
out in an adiabatic regime as the frequency of AF to is 
rather smallExLj. In this limit it was found for arbi- 
trary values of the noise strength and the amplitude of 
AF rj, the particle moves, in the direction of the slower 
rate of potential changaI3. However, in Ref. |ll| the pe- 
culiar reversal of the sign of v has been found as the AF 
frequency u) becomes relatively large. In this nonadia- 
batic limit and as the amplitude of AF rj increases, the 
dependence v{rj) displays oscillations. Moreover, these 
effects do not appear in the overdamped deterministic 
regime where the noise strength is zerocj. Notice here, 
that the directed transport of Brownian particle is more 
complex than the one in the deterministic case. In latter 
case, the non zero value of v appears as a simple conse- 
quence of two different values of a depinning force, and 
the directed transport is absent as the amplitude of AF 
is small. Thus, in order to explain various peculiarities 
of the overdamped Brownian particle motion we develop 
a theoretical analysis that goes beyond an adiabatic ap- 
proximation. 

In this Letter we present a consistent analytical ap- 
proach to the dynamics of an overdamped Brownian par- 
ticle in the presence of bothJ?P U (x) with the period a 
and harmonic AF, 77 cos(wt;)Ej. By making use of a par- 
ticular diagrammatic technique that is valid in the limit 
of a large time, the dependence of the mean velocity on 
the amplitude rj and frequency u> of AF will be calculated. 
We find that in the presence of AF the average velocity 
is determined by various relaxation processes inside the 
potential well, and the directed transport of Brownian 
particle occurs as the potential U (x) has a low symmc- 
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try. We also obtain the particular range of u> where the 
interplay between noise and AF leads to the reversal of 
the sign of v. 

The dynamics of an overdamped particle is described 
by the Langeven equationQ : 



ax(t) + U'(x) = r/ cos(ujt) + £(t) 



(1) 



where the white noise function £ (t) has a zero mean and 
the correlation function < £(i)£(£') > = 2aTS{t - t'). 
Here, a is the damping coefficient and T is the effective 
temperature describing the strength of the fluctuations. 
Next, we introduce the time-dependent probability idgn- 
sity P(x, t) that satisfies the Fokker-Plank equatiorc'B: 



particle return to the stable state after the multiple in- 
teractions with AF. Moreover, we are interested in a non 
oscillating part of P(x,t) only. It allows greatly to sim- 
plify the Eq. (^) just keeping most singular contributions 
with small "lambda" arguments of the Green functions 
G\{x,y). These terms can be presented in the diagram- 
matic form (see, Fig. 1). By summing these terms and 
calculating the intermediate integrals over A„ we obtain 
the expression for G\(x, xo) in the form: 

G x (x,x a ) = J dpe^-^GKp; x, x ) - ^ . (7) 
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Here, we introduce the periodic Green function 
Gl(p; x, xq) written in the mixed momentum -coordinate 
{[U'(x) - r]cos(uit)]P(x,t)} representation, and Z x (p) is determined by expression: 



with the initial condition 

P(x,0) = S(x-x a ) 



Using 



the 



Laplace 



(2) 



(3) 

transformation 



P\{x) = J Q dtP(x,t)e xt and the standard substi 

gfap — t/(»o) 

tution P\(x, Xo) = e ^ G x (x, Xq) we obtain the 
integral equation for the Green function G\{x, xq) 

d\x 2Ai 



G\{x,x ) = G° x (x,x a ) + / —77 



2iri 'u 2 + X 2 



/ dyG x (x,y)e *r — [e ™ Gx-Xt(y,x )} , 
a J dy 



(4) 



where G'Kx, Xq) is the Green function of the equation: 



Zx(p) 
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dyidy 2 P x (p; 2/2,2/1) 
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(8) 



T-^G x (x, x ) + G° x (x, x )(-Te^ -aX) = 

— ~aS(x — Xq) . (5) 

The mean value of particle velocity v is determined by 
the probability P(x, t) in the limit of large time: 



= lim 

t — >oo 



dxP{x, t) 



(x - x ) 



(6) 



In the absence of AF (77 =p-, 0) it immediately follows from 
Eq. (|) andthepropertjJiaG^xo) = G° x (x ,x) that 
the mean value of the velocity is zero for an arbitrary pe- 
riodic potential U{x). However, in the presence of AF the 
probability density P(x, t) is determined by the product 
of the Green functions with different arguments A. It 
leads to the breaking of the symmetry and, therefore, to 
the directed transport. 

To obtain the fluctuation induced transport in the 
presence of AF we need to know the long time behaviour 
of the probability distribution P(x,t) (see, the Eq. (||)). 
This limit corresponds to small values of A in the inte- 
gral equation M) . Physically it may be interpreted as the 
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FIG. 1. The most important diagrams and the diagram- 
matic presentation of Eq. (0) . Thin and thick solid lines cor- 
respond accordingly to G° x (~p\ x,y) and P x (x,y)\ dashed lines 
are due to the presence of AF. The cross presents an operator 
In the limit of large time the parameters p and A are 

small 



To calculate the function Z x (p) and correspondingly, 
to obtain the conditions when the directed transport oc- 
curs, we use the general properties of the Eq. (H). This 
equation can be mapped to a well known problem of an 
electron motion in a periodic potentiaKra. Thus, the re- 
laxation times spectrum T of the Eq. (^]) contains an 
infinite number of bands and is determined by the wave 
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vector p. The lowest band starts from the zero value, and 
in the region of small wave vectors p has a form: 
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2tt 



= Dp', p < 

Tn(P) a 



(9) 



where the diffusion coefficient D depends on the strength 
of the potential and temperature. In the limit of small 
fluctuations as the amplitude of U(x) is large and the 
temperature T is small, we obtain that the value of 

D oc e t = e~~ is also small. Note here, 

that the diffusion coefficient D does not depend on the 
symmetry of the potential U and is determined by the 
the narrow regions close to the top and bottom of poten- 
tial. In the absence of AF the particle displays a diffusive 
motion and it is completely determined by D. The lowest 
band pan be also mapped to a well known tight-binding 
modelD. The corresponding eigen function tjsp(x) has the 
form of a Bloch waveo, namely ip p (x) — e lpx r p (x), where 
the periodic function 



r p (x) ~ exp(- 



U[x) 
2T 



) + pu (x) 



(10) 



in the limit of small p. Here, the first correctk*-,uo (x) in 
the region of small p has been found in Refs.SO. Upper 
bands correspond to much smaller relaxation times r n 
but exactly these bands determine the directed transport. 

In the limit of large time the main contribution to 
the probability distribution P(x,t) results from a small 
momentum p. The presence of an odd component of 
Z\(jp) oc p leads to a nonzero value of v, and thus, 
by using @, @ and substituting (§) and (@) in (g) we 
arrive at the expression for v 
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(11) 



where A = J Q a dyi dy 2 e t i} e t 2 . Here, the Green 
function G\(x, xq) is a solution of the Eq. fl|) with peri- 
odic boundary conditions. 

The formula ( |Tl| ) is convenient for analysis of both gen- 
eral properties of Brownian particle motion and partic- 
ular limits. Thus, a general feature of an overdampcd 
Brownian particle transport is that the average velocity 
v is determined by the properties of potential U(x) on 
the whole range of its variation. Moreover, the directed 
transport of Brownian particle is an interplay of two ef- 
fects: the fluctuation induced particle escape from the 
potential well (the coefficient 1/A in the Eq. ([ll])) and 
the AF induced relaxation processes inside the potential 
well (the Green function Gi u (yx, j/2) in the Eq. (pT|)). 
However, if the potential U(x) displays a reflection sym- 
metry, namely U (x) is an even function of x, the average 
value of v vanishes due to a fact that -^Piu{yi,y%) is an 



odd function of y\ or y 2 , and one of the integrals in the 
Eq. (|Tl| ) is zero. 

If the potential U(x) has no reflection symmetry, the 
average velocity is not zero, and it increases as rj 2 in 
the limit of a small amplitude of AF rj. In this limit by 
making use of the eigen functions tjs n (x) of the Eq. (J5|) 
we get: 



< v > 



rj 2 2a \ - 

n 



' (^T„) 2 



dyi / dy 2 e ^ e ^ ^ n (y 2 )ip n (yi) . (12) 
Jo 

Thus, the average velocity v is small in both limits of 
small and large effective temperature T: 



v(T) oc 



e-^\ T < U. 



U„ 



T- 3 f(wn) , T » U max -U, 



max ^ in, in 
min 



(13) 



The particular function / depends on the ratio between 
the frequency lu and t^ 1 where t\ is the maximum re- 
laxation time of the particle in the potential U(x). The 
Eq. ( p"2| ) shows that the adiabatic regime is valid only in 
the limit ofw < 1 jr\ , The relaxation time T\ depends 
crucially on the effective temperature T. Thus, as the 
fluctuations are weak (T <C At/), T\ = where 

ujq is the frequepcy of small oscillations at the bottom 
of potential welllia. Moreover, n can be even larger in 
the limit of small fluctuations as thfi potential U{x) has 
a form of a double potential welli3'E£l. In the opposite 
regime of large fluctuations (T ^> At/) the relaxation 
time depends on T and decreases as 



n 



4tt 2 T 



(14) 



In the limit of a small frequency us the main contri- 
bution to the directed transport results from the first 
"relaxation time" band. It leads to a particular sign of 
v that does not change with the temperature. However, 
as the frequency increases (us 3> r^ 1 ) the upper "re- 
laxation time" bands dominate and the sign reversal can 
occur. As T increases the relaxation time t\ decreases 
and the adiabatic limit is recovered. 

This general scenario can be verified in the limit of 
large temperature T 3> At/ for a particular model of low 
symmetry potential U(x) = t/ (cos(^) + 7 sin(^)), 
where the parameter 7 is of order one. As the temper- 
ature is large the eigen functions ipn(x) of the Eq. (||) 
are plain waves, the main contribution to the v results 
from the first and second relaxation time bands, and we 
obtain: 



a z a 1 J 1 + (u!Ti) z 



1 



16 + (us Tl ) 2 



(15) 



where t\ is determined by the Eq. (|IJ) . The dependence 
of v(us) that displays a sign reversal, is presented in Fig. 
2. 
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In the discussion presented to this point we assume 
that the amplitude of AF r\ is small. As r\ increases the 
peculiar oscillations appear in the dependence of v{rj). 
The period of these oscillations A77 can be estimated in 
the same limit of large temperature T. By making use of 
the Eq. and a smallness of the potential U we find: 

A„(T) oc ( am i ' W » T U (16) 
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FIG. 2. The frequency dependence of the mean value of 
velocity v. the limit of large effective temperature. 

In conclusion, we have presented an analysis of an over- 
damped Brownian particle motion in the presence of pe- 
riodic potential and AF. We have shown that the directed 
transport in such a system is determined by two effects: 
the particle escape from a potential well and the relax- 
ation processes inside the potential well. It is at variance 
with both the standard diffusive motion in the absence of 
AF and tho.pDpJinear dynamic response of the system to 
a weak AFQO'Ej. We argue that the w dependent rever- 
sal of the sign of v found previously by numerical analysis 
is due to different strength and symmetry of relaxation 
processes. Thus, the observation of directed transport 
and especially its dependence on the frequency u> and 
amplitude 77 of AF, allows to investigate the various re- 
laxation processes of an overdamped Brownian particle. 
Moreover, the used method of taking into account the 
most singular diagrams where the momentum p and pa- 
rameter A are small, can be also useful in order to analyze 
the fluctuation induced transport in .the presence of more 
complex (periodic or random) AFE3l_3. 

I thank S. Flach, O. Yevtushenko, Y. Zolotaryuk and 
A. V. Ustinov for useful discussions. 
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